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Abstract 



The single-particle and collective dynamical properties of liquid lithium 
have been evaluated at several thermodynamic states near the triple point. 
This is performed within the framework of mode-coupling theory, using a self- 
consistent scheme which, starting from the known static structure of the liq- 
uid, allows the theoretical calculation of several dynamical properties. Special 
attention is devoted to several aspects of the single-particle dynamics, which 
are discused as a function of the thermodynamic state. The results are com- 
pared with those of Molecular Dynamics simulations and other theoretical 
approaches. 

PACS: 61.20.Gy; 61.20.Lc; 61.25.Mv 
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I. INTRODUCTION 



For the past twenty years or so, the dynamics of liquid metals has been a field of intense 
research, both theoretical and experimental (see, e.g., Ref. [I]), especially as regards the liquid 
alkali metals, all of which have been studied experimentally by means of inelastic neutron 
scattering (INS) or inelastic X-ray scattering (IXS) or by both techniques: LijHI Na,@~i 
K,0 RbfMl and Cs.0 Molecular Dynamics (MD) simulations have also stimulated this 
progress because of their ability to determine certain time correlation functions which are 
not accessible to experiment, and thereby they supplement the information obtained from 
the experiments. Most such MD studies have likewise been performed for alkali metals,EHli 
although other systems such as the liquid alkaline-earths!^ and liquid leadll have also been 
studied. 

On the theoretical side, this progress can be linked to the development of microscopic 
theories which provide a better understanding of the physical mechanisms underlying the 
dynamics of simple liquids.^ An important advance in this respect was the observation 
that the decay of several time-dependent properties can be explained by the interplay of 
two different dynamical processes.0IHIzl The first one, which gives rise to a rapid initial 
decay, comprises fast, uncorrelated, short-range interactions (collisional events) which can 
be broadly identified with "binary" collisions. The second process, which leads to a long- 
time tail, is connected to the non-linear coupling of the dynamical property of interest with 
slowly varying collective variables ("modes") such as density fluctuations, currents, etc., and 
is therefore referred to as a mode-coupling process. 

By introducing some simplifying approximations, Sjogren and coworkersiUnl first applied 
this theory to evaluate several collective and single-particle dynamical properties of liquid 
rubidium as well as the single-particle properties of liquid argon, at thermodynamic con- 
ditions near their respective triple points, obtaining results in qualitative agreement with 
the corresponding MD simulations. Balucani and coworkersEHl also applied a simplified 
mode-coupling theoretical approach to study the liquid alkali metals close to their triple 
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points, obtaining good qualitative predictions for several transport properties (self-diffusion 
coefficient and shear viscosity) as well as other single-particle dynamical properties (velocity 
autocorrelation function and its memory function, and mean square displacement). Based on 
these ideas, we developed a theoretical approach^ which allows a self-consistent calculation 
of all the above transport and single-particle dynamical properties; its application to liquid 
lithiumS and the liquid alkaline-earths^ near their triple points, has lead to theoretical 
results in fair agreement with both simulations and experiment. 

Since single-particle and collective dynamical properties are closely interwoven within the 
mode-coupling theory, the application of this formalism to any liquid system should imply 
the self-consistent solution of the coupled equations appearing in the theory. However, none 
of the above-mentioned theoretical calculations have been performed in this way. In fact, 
the usual practice is to obtain the input dynamical quantities needed for the evaluation of 
the mode-coupling expressions either from MD simulations or from some other theoretical 
approximation. In particular, this has often been the case for the intermediate scatter- 
ing function, F(k,t), which has been obtained either from MD simulations,0i3 or derived 
from some simple approximation such as the viscoelastic model or some simplification of 
it .SHIES 

We have recently presented a theoretical scheme^ for the self-consistent determination of 
F(k, t) within the mode-coupling theory. Its application to study the dynamical properties of 
liquid lithium at thermodynamic states near the triple point lead to good qualitative results 
in comparison with MD simulations. However, the self-intermediate scattering function, 
F s (k, t), was evaluated by means of the gaussian approximation, leading to an unsymmetric 
treatment of F{k,t) and F s (k,t), because mode coupling effects were included for the first 
but not for the latter. This imbalance is corrected in the present paper, in which we introduce 
a self-consistent framework which treats on an equal footing, and within the mode coupling 
theory, both F(k,t) and F s (k,t). This new scheme is then applied to study the dynamical 
properties of liquid lithium, covering a somewhat larger range of temperatures than in our 
previous study. 



The only input data required by the theory are the interatomic pair potential, which was 
obtained from the neutral pseudoatom (NPA)S method, and the liquid static structural 
functions, which were evaluated through the variational modified hypernetted chain approx- 
imation (VMHNC)&~0 theory of liquids. Comparison with MD simulations^! has shown 
that this combination leads to an accurate description of the equilibrium properties of liquid 
lithium close to the triple point. We stress that by combining the NPA method to obtain 
the interatomic pair potential, and the VMHNC to calculate the liquid static structure, we 
are able to attain the required input data, for the calculation of the dynamical properties, 
from the only knowledge of the atomic number of the system and its thermodynamic state. 

The paper is organized as follows. In section |I| we describe the theory used for the 
calculation of the dynamical properties of the system and we propose a self-consistent scheme 
for the evaluation of both single-particle and collective properties. In section |TXT| we present 
the results obtained when this theory is applied to liquid lithium at some thermodynamic 
states. Finally we sum up and discuss our results. 



II. THEORY 

This study is based on a combination of kinetic theory ideas with the Mori memory 
function approach. Within this framework, it is assumed that the memory functions 

of several time correlation functions are controlled by the interplay of two different dy- 
namical processes: one, with a rapid initial decay, due to fast, uncorrelated short-range 
interactions; the second (known as a mode-coupling process), with a long-time tail, due to 
the non-linear coupling of the specific time correlation function with slowly varying collective 
variables ("modes") such as density fluctuations, currents, etc. In the present theoretical 
framework, we focus on two basic dynamical variables, namely, the intermediate and the 
self-intermediate scattering functions which provide information about the collective and 
single-particle dynamical properties of the system, respectively. Moreover, their respec- 
tive time-Fourier transforms give the dynamic and self-dynamic structure factors which are 
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amenable of determination by means of both INS and IXS experiments. 



A. Collective Dynamics 

The collective dynamical properties are embodied in the dynamic structure factor, 
S(k,ui), which can be obtained as 

S(k, uj) = - Re F(k, z = -iu) , (1) 

7T 

where Re stands for the real part and F(k, z) is the Laplace transform of the intermediate 
scattering function, F(k,t), i.e., 



oo 

zt 



F{k,z) = / dte- zt F{k,t). (2) 



o 



Within the memory function formalism, F(k, z) can be expressed as 

n 2 (k) 



F(k,z) = S(k) 



z + 



(3) 



z + T(k, z 

where T(k,z) is the Laplace transform of the second-order memory function, T(k,t), and 
Q 2 (k) = k 2 / (3mS(k) where m is the mass of the particles, (3 is the inverse temperature 
times the Boltzmann constant and S(k) is the static structure factor of the liquid. Note that 
since we consider spherically symmetric potentials and homogenous systems, the dynamical 
magnitudes only depend on the modulus k =| k \. 



Now, the second-order memory function T(k,t) is decomposed as follows:oc3EZl 

r(M) = r B (M) + r M c(M) , (4) 

where the term r^/c, t), known as the binary part, has a fast decay, whereas the the mode- 
coupling contribution, Tmc^^), which aims to take into account repeated correlated col- 
lisions, starts as t 4 , reaches a maximum and then decays rather slowly. We now briefly 
describe both terms; for more details we refer the reader to Ref. |I|,|2"B|-p7j. 
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1. The binary- collision term 



Over very short times, the memory function is well described by T B {k,t) alone, both 
T(k,t) and T B (k,t) having the same initial value and curvature:0 

v(k, o) = r B (k, o) = ^- + nl + y d (*) - n 2 (k) , (5) 

where 



fi o = ^ I drg(r)V 2 <p(r) (6) 



is the squared Einstein frequency and 



l l d {k) = -JL( dre-^gir) (k ■ V)V(r) , (7) 



m 



with (p(r) and g(r) denoting respectively the interatomic pair potential and the pair distri- 
bution function of the liquid system with number density p, and k = k/k. 

The binary term includes all the contributions to T(k,t) to order t 2 . Since the detailed 
features of the "binary" dynamics of systems with continuous interatomic potentials are 
rather poorly known, we resort to a semi-phenomenological approximation that reproduces 
the correct short-time expansion. Therefore we writeS 

r B (k,t) = r B (k,o)e- t2 ^ k \ (8) 

where the relaxation time, Ti(k), can be determined from a short-time expansion of the 
formally exact expression of the binary term, and is related to the exact sixth frequency 
moment of S(k,u). In this way, after making the superposition approximation for the 
three-particle distribution function, one obtainsB 

r B (M) 



3k 2 



2Pm 
, 3p 



2k 2 

_ +3(^ + 2^) 



+ 



(ttt)^ / dfe- ik - r g(r)(k ■ V)V(r-) + 



6 




(9) 



where summation over repeated indices is implied (a, [3, 7 = x,y,z), and 



m 



P 



I 



(10) 



The relaxation time can thus be evaluated from the knowledge of the interatomic pair 
potential and its derivatives together with the static structural functions of the liquid system. 



The inclusion of a slowly decaying time tail is essential in order to achieve, at least, a 
qualitative description of T(k, t). A detailed treatment of this term requires consideration of 
several modes (density-density coupling, density- longitudinal current coupling and density- 
transversal current coupling). However, for thermodynamic conditions near the triple point, 
the most important contribution is the density-density term, which can be written as0 



where c(k) is the direct correlation function and F B (k, t) denotes the binary part of the 
intermediate scattering function, F(k,t). Following Sjogren, HI we approximate the ratio be- 
tween F(k, t) and its binary part by the ratio between their corresponding self counterparts, 
i.e., 



2. The mode- coupling component 




x 'k-k'c(k') + k-(k-k')c(\ k-k'\) 



x F(\ k-k' \,t)F(k',t)-F B (\ k-k' \,t)F B (k',t) 



(11) 



F B (k,t) 



F s , B (k,t) 
F s (k,t) 



F(k,t), 



(12) 
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where F s (k, t) is the self- intermediate scattering function and F S) B{k, t) stands for its binary 
part, which is now approximated by the free-particle expression 

F StB (k,t) = F (k,t) = ex V [-^k 2 t 2 } . (13) 

Once F s (k,t) has been specified, the self-consistent evaluation of the above formalism 
will yield to F(k,t), and, therefore to all the collective dynamical properties discussed so 
far. 

3. Other models 

In this work, we will compare the predictions of the above self-consistent approach with 
those of a much simpler and widely used model for the second-order memory function, 
T(k,t); this is the so-called viscoelastic model. It approximates T(k,t) by an exponentially 
decaying function with a single relaxation time, which is usually fitted so that the predicted 
value of S(k,u = 0) coincides for k — > oo with the exact free-particle result.0@ 

B. Single-particle dynamics 

The self-intermediate scattering function, F s (k,t), probes the single-particle dynamics 
over different length scales, ranging from the hydrodynamic limit (k — > 0) to the free-particle 
limit (k — > oo). Its frequency spectrum is the self-dynamic structure factor, S s (k,u) = 
(l/7r) Re F s (k,z = —ioo), where F s (k,z) stands for the Laplace transform, which according 
to the memory function formalism, can be expressed as 



FJk.z) = : + h'Ak.z) _1 = z+ 



(14) 



z + T s (k, z) 

where ^(k) = k 2 /(j3m) and K s (k,z) and T s (k,z) are respectively the Laplace transforms 
of the first- and second-order memory functions of F s (k,t). Note that the velocity autocor- 
relation function (VACF) of a tagged particle in the fluid, i.e., | (vx(t)vi(0)) , is also given 
by the limit linu^o K s (k, t)/k 2 . As before, T s (k,t) is now decomposed!!!! 



r s (k, t) = r S}B (k, t) + T sMC (k, t) , 

into a binary contribution, T StB (k,t) and a mode-coupling term, T Sj uc(k,t). 



(15) 



1. The binary- collision term 

The binary term shows a fast decay, T s (k,t) and T Sj B{k,t) having the same initial value 
and curvature :0 

2k 2 

T s (k,0)=T s , B (k,0) = — + n 2 . (16) 

pm 

Its time dependence can be adequately described by a semi-phenomenological expression, 

r SiB (M) = r s ,£(M)e^ 2/Tfl2(fc) , (17) 

where T s (k) is a relaxation time, determined from a short time expansion of the formally 
exact expression of the binary term, which can be related to the sixth frequency moment of 
S s (k,u). In this way, after making the superposition approximation for the three-particle 
distribution function, one obtains^ 



r fliB (fc,0) 3A; 



2 



r 2 (k) 2(3m 



2 



2/c o 

— + 3n 2 

pm 



2 



+ -3 (is) 



r 2 



where 



t 2 3m 2 



df[V a V^{r)} g(r) [V Q V /3 ^(r)] 



1 r dk 1 

+ 67 / (2^ 7 ^' } ^ ~ > (19) 

with summation over repeated indices implied. 

2. The mode- coupling component 

This term describes the intermediate and long time behaviour, involving repeated corre- 
lated collisions. As for its collective counterpart, we simplify the full multi-mode expression 
by considering only the coupling to the density fluctuations^ 



T sM c{k,t) = 

x 



dk' 



(k-k') 2 c 2 (k" 



f3m J (2tt) 3 

F s (\ k-k' \,t)F(k',t)-F s>B (\ k-k' \,t)F B (k',t) 



(20) 



Note that the term involving the product F s ^Fb, which makes T s ^c{k,t) very small at 
short times, decays very fast. Therefore the possible errors introduced by making specific 
approximations for these binary terms become negligible in the region of intermediate and 
long times, where the dominant contribution to T sMC (k,t) is the product F S F. 



3. Other models 

A viscoelastic model has also been proposed for the F s (k,t); it assumes for T s (k,t) an 
exponentially decaying function with a single relaxation time, which can in principle be fixed 
using several prescriptions that have been proposed. In this work we will use the expression 
obtained bu requiring that S s (k,u = 0) coincides for k — > oo with the exact free-particle 
result .0 

A different, more accurate and widely used, model for F s (k,t) is the so-called gaussian 
approximation (GA): 

1 



F s (k,t) = exp 



-- k 2 Sr 2 (t) 
6 V ; 



(21) 



where 5r 2 (t) = (\fi(t) — n(0)| 2 ) stands for the mean squared displacement of a tagged 
particle in the fluid. This approximation produces correct results in the limits of both small 
and large wavevectors, and for all wavevectors at short times. Moreover, 5r 2 (t) is related to 
the normalized VACF, Zit) = (v^v^O)) / (v 2 ) , by 

5r 2 {t) = ^-fdr{t-r)Z{r) . (22) 

Lj lib J 
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C. Self-consistent procedure 



The calculations begin with some estimation for both F(k,t) and F s (k,t) (e.g., those 
estimates provided by the corresponding viscoelastic models). Using the known values of 
^s,B{k,t) and Eq. (^Dj) and (|T2|), a total memory function T s (k, t) is obtained which, when 
taken to Eq. (|HD, gives a new estimate for F s (k,t). This is now taken to Eq. ( |TTD which, 
along with the known values of Ts{k, t), produces an estimate for the total memory function 
T(k,t), which is taken to Eq. (^) for a new determination of F(k,t). Now, with the new 
estimates for both F{k,t) and F s (k,t) the whole procedure is repeated. 

The previous computational loop is iterated until self-consistency is achieved between 
the initial and final F(k, t) and F s (k, t). The practical application of this scheme has shown 
that reaching self-consistency requires about ten iterations. 



D. Transport properties 



Once the self-consistency has been achieved, the normalized VACF is obtained as 



Z(t) = ((3m) Ihn^M (23) 



whereas its associated transport coefficient, namely the self-diffusion coeficient, D, is given 
by 

1 r°° 

D = — dtZ(t). (24) 
pm Jo 

Another interesting transport property is the shear viscosity coefficient, r], which can 
be obtained as the time integral of the stress autocorrelation function (SACF), r](t), which 
stands for the time autocorrelation function of the non-diagonal elements of the stress tensor. 
Moreover, r](t) can be decomposed into three contributions, a purely kinetic term, 
a purely potential term, rj pp (t), and a crossed term, rjk p (t). However, for the liquid range 
close to the triple point, the contributions to rj coming from the first and last terms are 
negligible,il and therefore in the present calculations we assume r/(t) = f] pp (t). This function 
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can in turn be split into its binary and mode-coupling components, r](t) = T)B(t) + f]Mc{t)- 
Again, the binary part is described by means of a gaussian ansatz, i.e., 



(25) 



where the rigidity modulus, G p , and the initial time decay, can both be computed from the 
interatomic potential and the static structural functions of the system.0 The superposition 
approximation for the three-particle distribution function is also used in the evaluation of 
t v . For the mode-coupling component, r]Mc(t), we only consider the coupling to density 
fluctuations, given by,0BS 



VMc(t) 

1 

60/?7T 2 



dkk 4 



S'{k) 
S 2 (k) 



F 2 (k,t) - F 2 B (k,t) 



(26) 



where S'(k) is the derivative of the static structure factor with respect to k. This mode- 
coupling integral is evaluated using the F(k,t) and Fs(k,t) previously obtained within the 
self-consistent scheme described above. 



III. RESULTS 

We have applied the preceeding theoretical formalism to study the dynamical properties 
of liquid 7 Li at several thermodynamic states near the triple point (see Table I). The input 
data required for the calculation of the dynamical properties are both the interatomic pair 
potential and its derivatives, which were evaluated by the the npaH method, as well as the 
liquid static structural properties, evaluated by the VMHNC§@ theory of liquids. 

First, we have evaluated the relaxation times appearing in the binary part of the second 
order memory functions of the intermediate and self-intermediate scattering functions, Ti(k) 
and r s (k), respectively, using Eq. (^) and (|TB|). The obtained results are plotted in Fig. [I], 
which shows that for each temperature, both relaxation times have approximately the same 
magnitude, with ri{k) oscillating around r s {k). 
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A. Intermediate scattering function 



Figures |2| and |3] show, for T=470 K and 725 K, at two k values, the first and second-order 
memory functions of F(k,t), as obtained from the MD simulations, and the corresponding 
theoretical ones derived within the present approach. The MD results for the second-order 
memory function, Fmd^^); show a rapid initial decay followed by a long-time tail which 
becomes smaller when the wavevector or the temperature increases, and exhibits an oscil- 
latory behaviour for k > k p /2. The MD results for the first-order memory function show 
a negative minimum whose absolute value decreases as k or T increases. The theoretical 
T(k,t) obtained by the self-consistent procedure qualitatively reproduce the corresponding 
MD results, especially the short time behaviour, which is dominated by the binary compo- 
nent. For larger t the overall amplitude of the decaying tail of T(k,t) is, in general, well 
described, but some discrepancies with the MD results appear for the amplitudes of the os- 
cillations. We believe that its improvements will require the inclusion of other terms in the 



expression for TMc(k,t) (see Eq. (|TTD), in particular those related to the density-currents 
couplings. The viscoelastic model,i! whose results, r visc (fc,t), are also shown in Figs. ^| and 
|3], is clearly unable to describe the features exhibited by TyiD(k,t). 

The corresponding MD results for the intermediate scattering functions, -Fmd (.k, t) , are 
shown in Figs. ^ and [5|. Below k pa kp/2, FuD{k,t) exhibit an oscillatory behaviour, with 
the amplitude of the oscillations being stronger for the smaller /c-values. These features 
are predicted quite well by the self-consistent theoretical F(k,t), especially as regards the 
amplitude of the decaying tail, which is well reproduced. The main discrepancy is the 
underestimation of the amplitude of the oscillations for small k. These results, however, 
represent an important improvement over those obtained from the viscoelastic model,!! 
F V i SC (k,t), which basically oscillate around zero. Moreover, for small fc-values, i.e., k ~ 0.25 
A -1 , the oscillations in F v i sc (k, t) are overdamped, whereas for the intermediate k- values, i.e., 
k pa 1 A -1 and k ~ 1.7 A -1 , they are too strong. Only for A;-values around k p the viscoelastic 
results show a good agreement with the MD ones. This agreement can be understood as a 
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compensation between the short-time and intermediate-time deficiencies of T v i SC (k,t). On 
the other hand, note that the viscoelastic model imposes the exact initial values of F(k,t), 
and its second and fourth derivatives; this explains the good agreement obtained at short 
times. 

Finally, we point out that the theoretical results for F(k, t) and its corresponding memory 
functions, are rather similar to the those obtained using the GA model for F s (&,£). This 
is not unexpected, because the role of F s (k, t) in the determination of T(k, t) is restricted to 



the evaluation of FB(k,t) appearing in Tuc(k,t) in Eq. (11). However, the part involving 
the product of the Fs(k, t)'s just makes Tuc{k,t) very small at short times, where, in any 
case, the dominant contribution is the binary term, r^(fc,t). Therefore the possible errors 
in F s (k,t) due to the GA are of little account. 



B. Self intermediate scattering function 

Figures ^| and show, for T=470 and 725 K, the MD and theoretical results obtained 
for the first and second-order memory functions of F s (k,t). 

The MD results for the second-order memory function, ^ s ,MD{k, t), are qualitatively very 
similar to those previously obtained for Fmd(^> t): it exhibits a rapid decaying part at short 
times, followed by a long-time tail which can take negative values for the greater fc-values. 
However, the long-time tail of T Sj ui){k,t) oscillates for all fc-values, not just for k > k p /2, 
as in the case of TuD{k,t). Also, we note that the role of the long-time tail becomes less 
relevant when the temperature is increased. The first-order memory function has a negative 
minimum for all the fc-values, followed by an oscillating tail. In fact, similar qualitative 
features were also obtained by Shimojo et a/.0 in their MD study for liquid Na near the 
triple point. The self-consistently calculated T s (k,t) shows a good overall agreement with 
the MD results. The short-time behaviour, which is dominated by the binary component, 
is well reproduced, since the present theoretical approach imposes the exact initial values of 
both T s (k,t) and its second derivative. The intermediate and long-time behaviour, which is 
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controlled by F St uc(k,t), is qualitatively well described, especially the oscillations of the tail 
for small fc-values. However, increasing the temperature (Fig. [7]) leads to an overestimation 
of F s (k,t) for intermediate and long times. This fact may signal the need to incorporate 
other coupling terms in the expression of r Sj Mc(^ ; t). For comparison, we have also included 
in Figs. |B| and [Tj the results of the viscoelastic model0 for the second-order memory function, 
r Si visc(&, t). It is observed that this function provides a rather poor description. In fact, for 
k < k p r S)Visc (fc, t) underestimates the MD results for all times, whereas for larger /c-values 
revise is too small for short times, and too large at the intermediate times. A further 
comparison is also performed with the GA model, as given by Eq. fl2lD , using the 5r 2 (t) 
obtained from the MD simulations. We have calculated the corresponding second-order 
memory functions, T sg (k,t), which are also shown in Figs. |6] and [j]. It is observed that 
r s>g (k,t) gives a reasonable account of the corresponding MD results, especially for short 
and intermediate times, whereas for longer times the model underestimates the MD results. 

The results obtained for the self-intermediate scattering function are shown in Figs. |8| 
and ^| for several k- values and T=470 and 725 K. The MD results, F SjMD (k,t), decreases 
monotonically with time for all the fc-values, and this behaviour is rather well described 
by the present theoretical formalism, which leads to F s (k,t) which closely follow the corre- 
sponding MD results. On the other hand, the viscoelastic model leads to F s y lsc (k,t) which 
underestimate the MD results seriously, especially for the smaller fc-values. This can be 
explained because the viscoelastic model incorporates the exact initial values of F s (k,t), 
and its second and fourth derivatives, but the relaxation time is fitted to produce the exact 
area of F s (k,t), i.e. S s (k,u = 0), for large k, whereas for low k the correct area is given 
in terms of the diffusion coefficient, which does not appear in the parametrization used. 
It can also be noticed that F SjV - lsc (k,t) does not vary monotonously with time for large k, 
showing oscillations which do not appear in the MD or in the self-consistent results. The 
GA model leads to F S:g (k, t) which compare favorably with the MD results, although for the 
intermediate fc-values predicts a quicker decrease. Taking z = in Eq. (|TJ|) shows that this 
is a consequence of the previously mentioned underestimation of the corresponding T s>g (k, t) 
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at long times. 

By Fourier transforming F s (k, t) we obtain the spectrum S s (k, u) which, for all /c-values, 
exhibits a monotonic decay with frequency, from a peak value at u — 0. In fact, the 
relevant features embodied in S s (k, to) are conveniently expressed in terms of the peak value 
S s (k,u = 0), and the half-width at half maximum, uoi/ 2 {k). These magnitudes are usually 
reported normalized with respect to the values of the diffusive (k — » ) limit, introducing 
the dimensionless quantities = nDk 2 S s (k,uj = 0) and A(/c) = uJi/ 2 {k)/Dk 2 . The 

magnitude uj\/2{k)/k 2 can be also interpreted as an effective fc-dependent diffusion coefficient 
D(k). For a liquid near the triple point, A(/c) usually exhibits an oscillatory behaviour 
whereas, in a dense gas it decreases monotonically from unity at k = to the 1/k behaviour 
at large k. The results obtained for A(k) in this work are shown in Fig. [IIJ The MD 
values, Amd(^), show that for all temperatures, the diffusive limit is reached from below, 
with a minimum at around k ~ k p , followed by a maximum and by a gradual transition, 
for greater k- values, to the free-particle limit. This oscillating behaviour of Amd(^) for 
small and intermediate fc-values has already been studied by several authors00il'@~§! and 
has been attributed to the coupling of the single-particle motion to other modes in the 
system; in terms of the present theoretical formalism this effect would be described by the 
F St Mc(k,t) term in Eq. (|15|). We note that similar features to those obtained in this paper 
for Amd(^); were already obtained by Torcini et a/.0 in their MD study of liquid lithium 
using the interatomic pair potentials proposed by Price et a/.0 

The self-consistent theoretical results obtained for A(k) exhibit a qualitative agreement 
with the corresponding MD ones. In particular the positions of the minimum and maximum 
of Amd(^), are successfully reproduced, although they are somewhat overeshot. These 
features are closely related to the term r SjMC (/c, t), which gives the intermediate and long 
time contributions of the second order memory function. This is shown by the fact that 
when only the binary term, T S) B{k,t), is considered in Eq. (p~5|) , the resulting values for 
the diffusion coefficient are poor, and the corresponding A#(A;) initially increases with k, 
contrary to the behaviour of the MD simulations Therefore, the inclusion of a tail in 
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T s (k,t) seems important in order to obtain the oscillating structure of Amd(&)- 



The self-consistently calculated values of are plotted in Fig. [11], where the compari- 
son with the MD results shows that the present formalism provides a rather good description 
of this magnitude. The shape of is more sensitive to changes in temperature, with the 
diffusive limit reached from below for T = 470 K and from above for T=725 K and 843 K, 
with the transition being located somewhere around T= 574 K. 

We point out that the present theoretical results for A(k) and £(&) are, to our knowledge, 
the first that have been obtained within a memory function/mode-coupling formalism, with 
no recourse to parameters, or inclusion of MD data, at any stage of the calculation. 



For comparison, we have also included in Figs, pfl and 11 the results obtained using 



the GA model of Eq. ([21]) for F s (k,t), using the MD values for Sr 2 (t). This approximation 
is exact for both small and large fc-values, whereas for intermediate fc-values the decay of 
F s (k,t) is too fast, producing too large a width and a smaller initial value of the spectrum 



S s (k, lu). This is clearly observed in Figs. |10| and |TT[ In fact, the associated A g (k) is unable 
to follow the fc-dependent behaviour exhibited by the corresponding MD results, especially 
the minimum appearing at k « k p . Note that the largest discrepancies appear for the inter- 
mediate fc-values, for which the spatial correlations are stronger, and that the discrepancies 
become smaller as the density is reduced. On the other hand, the GA results ontained for 
^g(k) qualitatively reproduce the MD results, although underestimating them. Moreover, 
as the temperature increases the agreement is improved. It is interesting to mention that 
although the A g (k) does not exhibit the oscillatory behaviour of Amd(^), however the cor- 
responding second-order memory function, T St9 (k, t) does have a tail which fairly follows the 
intermediate and long time behaviour displayed by T Sj uT)(k,t), as is observed in Figs. || 
and 0. Therefore, the existence of a tail in the second order memory function of F s (k,t), 
although necessary, does not automatically imply the appareance of an oscillatory behaviour 
of A(jfe). 

To gain further insight into the role played by the tail of T s (k,t) in the oscillatory 
behaviour of A(k), we have also evaluated F s (k,t) using an extension of the GA model 
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which includes some non-gaussian corrections, obtained by means of a cumulant expansion.El 
Restricting up to the first non-gaussian term, which provides the dominant corrections to 
the gaussian result, we have the following expression 



F s {k,t) = 



- k 2 5r 2 (t) 2 
L 6 W 



exp 



— - k 2 5r 2 (t) 
6 V ; 



(27) 



where a^t) is the first non-gaussian coefficient: 

3 5r\t) 



a 2 (t) 



(28) 



5 [5r 2 (t)} 2 

and 5r 4 (t) = (\ fi(t) — r*i(0) | 4 ). By using as input data the 5r 2 (t) and 8r 4 {t) obtained 
from the MD simulations, we have evaluated the corresponding F s (k,t), according to Eq. 
(fZTj). The associatted A n _ g (/c) and S n _ g (fc) accurately reproduce the corresponding MD 



results, as shown in Figs. |10| and [TTj. Moreover, the corresponding second-order memory 
functions, T Syn _ g (k, t), reproduce the MD results quite accurately. This is shown in Fig. |12 
for T = 470 K, where we have focused on the intermediate and long time behaviour, and 
the comparison is also performed with the MD simulations, the GA model and the present 
theoretical framework. Note that whereas the GA model always underestimates the MD 
results, the opposite happens with the present theoretical framework. The underestimation 
induced by the GA model is more marked for those /c-values around the main peak of S(k) 
which is precisely the region where Amd(^) exhibits a minimum. In order to explain the 
conection between this behaviour and the shape of the corresponding A(A;) = D(k)/D, we 
must note that D{k) = D(k, z = 0), where 

n 2 (k) 



k 2 D(k,z) = K s (k,z) 



(29) 



z + T s (k, z) 

where D(k, z) stands for a generalized diffussion coefficient.0 When this equation is taken 
at z = we obtain that D(k) T s (k,z = 0) = (l//3m). Therefore, according to Fig. [L2|, 



the GA model leads to a second order memory function, r ss (/c,t) which underestimates 
the corresponding MD results for all /c-values, and by the previous relationship it leads to 



greater estimates of the corresponding D(k), when compared with the corresponding MD 
results. By contrast, the opposite behaviour is exhibited by the present theoretical formalism 
as it overestimates the corresponding T s (k,t) and therefore leads to smaller values of the 
corresponding D{k). 



C. Transport properties 

The normalized velocity autocorrelation functions obtained from the self-consistently 
calculated results, using Eq. (p3|) , are shown in Fig. [13], where they a compared with the 
corresponding MD results. The initial decay of Z(t) is very well reproduced, because its 
initial value and both the second and fourth derivatives are implicitly imposed through Eq. 
([H]). The positions of the maxima and minima of Z(t) are also well predicted, although the 
amplitude of the oscillations is underestimated. The corresponding self-diffusion coefficients, 
D, are shown in Table [TT] along with those obtained by the MD simulations. For all the 
temperatures considered the MD results agree rather well with the experimental INS0@ and 
tracer data@ (for T = 843 K the experimental data have been extrapolated slightly outside 
the range suggested in Refs. |] and §). This good agreement supports the adequacy of the 
NPA-derived interatomic pair potentials to describe liquid lithium in this temperature range. 
The theoretical results show good agreement with the MD values for T=470 K, whereas as 
the temperature is increased the MD results are underestimated. In fact, according to the 
relation between D and the VACF, see Eq. (f24|), the smaller values obtained for D are a 
consequence of the above mentioned underestimation for the amplitude of the oscillations in 
Z(t), especially that of the first maximum. This is more evident at T=725 K, whereas for T 
= 470 K there is a cancellation between the first minimum and the first maximum leading 
to a theoretical D very close to the MD result. 

As regards the shear viscosity, fig. |H] shows, for T=470 and 725 K, the MD results 
obtained for the SACF along with its three contributions. Note that rjk p (t) and r)kk(t) parts 
of the SACF are very small (less than 10% of the potential-potential part) which justifies 
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their being ignored in the theoretical calculations of the SACF. The figure also shows the 
theoretical rj(t) along with its binary and mode-coupling components. Note that whereas 
for the lower temperature the inclusion of the mode-coupling component is essential for 
achieving a good overall agreement with the MD data, when the temperature is increased 
the role of the binary part becomes more dominant; in fact, for 725 K the binary part alone 
accounts for most of the MD results. 

Our present theoretical results slightly overestimate the short time (t ~ 0.1 ps) behaviour 
of T](t) and this is mainly due to the binary component; more explicitly, it comes from the 
overestimation of the values for t v . This limitation can be traced back to the use of the 
superposition approximation in the evaluation of the three-body term appearing in t v , which 
has a rather important weight in the final value.@ 

On the other hand, the behaviour for t > 0.1, which is completely determined by the 
mode-coupling component, is rather well described for both temperatures. The values ob- 
tained for the shear viscosity coefficient, rj, are presented in Table [TTI|. Although the the- 
oretical values slightly overestimate the MD ones (because of the short time behaviour of 
T](t)), both the theoretical and the MD results agree well with the experimental values.^ 



IV. CONCLUSIONS 

In this work we have evaluated several dynamical properties of liquid lithium at thermo- 
dynamic states close to the triple point. The calculations were performed by a self-consistent 
theoretical framework which, by incorporating mode-coupling concepts, allows the evalua- 
tion of both single-particle properties, as represented by the self-intermediate scattering 
function, its memory functions, the velocity autocorrelation function and the self-diffusion 
coefficient, as well as collective properties such as the intermediate scattering function, its 
memory functions, the autocorrelation function of the non-diagonal elements of the stress 
tensor, and the shear viscosity coefficient. Its application to liquid lithium has led to rea- 
sonable results when compared with both the corresponding MD results and the available 
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experimental data. The agreement is particularly satisfactory near the melting point, dete- 
riorating somewhat at higher temperatures. 

Within the present theoretical formalism, memory functions play a key role, specifically 
the second order memory function of both the intermediate scattering function and its self 
counterpart, as defined in Eq. (|3p and (0). All the relaxation mechanisms controlling 
the collective and single-particle dynamics, are introduced at the level of the corresponding 
second order memory functions which are decomposed into their binary and mode-coupling 
contributions. We have found that the gaussian ansatz adopted for the binary term pro- 
vides a reasonable description of the corresponding MD results for short times. As for the 
other contribution, namely the mode coupling part, for simplicity we have only considered 
the density-density coupling term. Although other couplings could be included, this term 
provides the dominant mode coupling contribution for thermodynamic states close to the 
triple point, which explains why both T(k,t) and T s (k,t) are better described for the lower 
temperatures. 

The self-consistent theoretical calculations provide results for T(k,t) and F(k,t) in qual- 
itative agreement with the MD results, and much more accurate than other theoretical 
approaches, as for example the widely used viscoelastic model. 

Better results have been achieved for T s (k,t) and the corresponding F s (k,t). Within 
this context, we emphasize the good description of the wavevector dependence of both the 
peak height and the half-width at half maximum of the S s (k,u), as represented by 
and A(k) respectively. Since these funcitons constitute a stringent test of any theoretical 
model,S we conclude that the main physical effects behind F s (k,t) seem to be included in 
the present theoretical framework, specifically, in the expression adopted for T s (k,t). We 
stress that, to our knowledge, this is the first theoretical study where the behaviour of A(k) 
has been qualitatively reproduced from just its memory function/mode-coupling formalism, 
without resorting to parameters, fitting to an assumed shape or including magnitudes from 
MD simulations in the evaluation of F s (k,t). 

The calculations carried out for the T s (k,t) have confirmed that the existence of a tail 
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is necessary to account for the oscillatory behaviour of A (A;); however, the magnitude of 
the tail, understood in the sense of its time integral, has a most important influence on the 
amplitude of the oscillations of A(fc). Moreover, the results obtained for the GA model with 
and without non-gaussian corrections, point towards the existence of a minimum magnitude 
of the tail in order to induce an oscillatory behaviour on A(k). 

The improvements achieved in the description of F s (k,t) and related magnitudes, when 
compared with those predicted by the GA model,0 are not fully reflected in the obtained 
values for the self-diffusion coefficients. However, this is not surprising because they are 
defined as a time integral of a correlation function which gives a measure of its time average 
but provides very scarce information on the dynamics of the system. 

We end up by signaling some limitations of the formalism presented here. First, its 
density/temperature range of applicability lies within the region where the relevant slow 
relaxation channel is provided by the coupling to density fluctuations, and this ceases to 
be valid for densities smaller than those typical of the melting point. For these densities, 
coupling to other modes, like longitudinal and/or transverse currents, becomes increasingly 
important and we believe that this is the main reason for the small deviations observed in 
the VACF total memory function at the two higher temperatures studied. In fact, further 
improvements in the description of the F s (k,t) and F(k,t), at all temperatures, would also 
require the inclusion of other modes in the corresponding second-order memory function. 
Moreover, although the previous remarks concern the mode-coupling contribution of the 
second order memory function, attention should also be drawn to the binary term. In fact, 
this term is poorly known and its role becomes increasingly dominant with large k and/or 
at increasing temperatures when the mode-coupling contribution decreases. A first task 
would be to study the influence of the superposition approximation for the three particle 
distribution function, on the obtained results for the relaxation times Ti(k) and T s (k). 

Further work is currently performed in that direction, and the results will be reported 
in due time. 
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FIGURES 

FIG. 1. Relaxation times ri{k) and T s (k), for liquid lithium. Continuous line and dotted line: 
ri(k) and T s (k) for T= 470 K. Dashed line and dash-dotted line: Ti{k) and r s (k) for T= 725 K. 

FIG. 2. Normalized second-order memory function, T(k, t), of the intermediate scattering func- 
tion, F(k, t), at two k- values, for liquid lithium at T = 470 K. Open circles: MD results. Continuous 
line: present theory. Dash-dotted line: binary part, Ts(k,t). Dotted line: mode-coupling part, 
Tyic(k,t). Dashed line: viscoelastic model. The inset shows the normalized first-order memory 
function as obtained by MD (open circles), the viscoelastic model (dashed line) and the present 
theory (continuous line). 

FIG. 3. Same as the previous figure but for T = 725 K. 

FIG. 4. Normalized intermediate scattering functions, F(k,t), at several A:- values, for liquid 
lithium at T = 470 K. Open circles: MD results. Continuous line: present theory. Dashed line: 
viscoelastic model. 

FIG. 5. Same as the previous figure but for T = 725 K. 

FIG. 6. Normalized second-order memory function, T s (k,t), of the self intermediate scattering 
functions, F s (k,t), at several /c-values, for liquid lithium at T = 470 K. Open circles: MD results. 
Continuous line: present theory. Dotted line: mode-coupling part, r Si Mc(k,t). Dashed line: 
viscoelastic model. Dash-dotted line: GA model for F s (k,t). The inset shows the normalized 
first-order memory function as obtained by MD (open circles), the viscoelastic model (dashed line) 
and the present theory (continuous line). 

FIG. 7. Same as the previous figure but for T = 725 K. 
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FIG. 8. Self intermediate scattering functions, F s (k, t), at several /c-values, for liquid lithium at 
T = 470 K. Open circles: MD results. Continuous line: present theory. Dashed line: viscoelastic 
model. Dash-dot line: GA model 

FIG. 9. Same as the previous figure but for T = 725 K. 

FIG. 10. Normalized half width of S s (k,uj), relative to its value at the hydrodynamic limit, for 
liquid lithium at four temperatures. Open circles: MD results. Continuous line: present theory. 
Dashed line: GA model for F s (k,t) . Dotted line: GA model with non-gaussian corrections for 
F s (k,t). Long dashed line: Free particle limit. Dash-dotted line: present theory with only the 
binary term in Eq. ( |T5| ) 

FIG. 11. Normalized peak value S s (k,uj = 0), relative to its value at the hydrodynamic limit, 
for liquid lithium at four temperatures. Open circles: MD results. Continuous line: present theory. 
Dashed line: GA model for the F s (k,t) . Dotted line: GA model with non-gaussian corrections 
for F s (k,t). Long dashed line: Free particle limit. Dash-dotted line: present theory with only the 
binary term in Eq. ( |l5| ) 

FIG. 12. Second-order memory function, T s (k,t), of the self intermediate scattering functions, 
F s (k,t), at several /c-values, for liquid lithium at T = 470 K. Open circles: MD results. Continuous 
line: present theory. Dash-dotted line: GA model for F s (k,t). Dotted line: GA model with 
non-gaussian corrections. 

FIG. 13. Normalized velocity autocorrelation functions of liquid Li at three temperatures. Open 
circles: MD results. Continuous lines: self-consistent calculations. 

FIG. 14. Normalized potential part of the stress autocorrelation function, rj(t), for liquid lithium 
at T=470 and 725 K. Open circles: MD results. Continuous line: theoretical results. Dash-dotted 
line: binary part. Dotted line: mode-coupling component. The inset shows the MD results for rj(t) 
(open circles), r] pp (t) (continuous line), 10 x rjk p (t) (dotted line) and 10 x r]kk(t) (dashed line). 
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TABLES 



TABLE I. Thermodynamic states studied in this work. 



T (K) 470 574 725 843 

p (A" 3 ) 0.0445 0.0438 0.0420 0.0416 



TABLE II. Self-diffusion coefficient (in A 2 /ps units), of liquid lithium at the thermodynamic 
states studied in this work. D t ^, and Dmd are the theoretical, viscoelastic and Molecular Dynamics 
results obtained in this work. 



T(K) 



470 



574 



725 



843 



Ah 



0.65 

0.69 
0.64±0.1 a 
0.69±0.12 b 
0.67±0.06 c 



0.99 

1.11 
1.08±0.15 a 
1.19±0.20 b 
1.16L0.09 c 



1.59 

1.94 
1.76±0.25 a 
1.99±0.30 b 
1.96±0.2 d 



2.02 

2.47 
2.28±0.30 a 
2.62±0.30 b 



a Ref. 
b Ref. 



c Ref. gS 
d Ref. M 
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TABLE III. Shear viscosity (in GPa ps) of liquid lithium at the thermodynamic states studied 
in this work. Tfth and t]md are the theoretical and Molecular Dynamics results obtained in this 
work. 



T (K) 470 574 725 843 

Tfth 0.59 0.46 0.36 0.29 

?7md 0.55 0.42 0.33 0.28 

i] cxp 0.57±0.03 a 0.45±0.03 a 0.35±0.03 a 0.30±0.03 a 

a Ref. M 
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